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Re-examining Archie’s law: Conductance description by tortuosity and constriction 
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In this article we investigate the electrical conductance of an insulating porous medium (e.g., a 
sedimentary rock) filled with an electrolyte (e.g., brine), usually described using the Archie cemen¬ 
tation exponent. We show how the electrical conductance depends on changes in the drift velocity 
and the length of the electric field lines, in addition to the porosity and the conductance of the 
electrolyte. We characterized the length of the electric field lines by a tortuosity and the changes 
in drift velocity by a constriction factor. Both the tortuosity and the constriction factor are de¬ 
scriptors of the pore microstructure. We define a conductance reduction factor to measure the local 
contributions of the pore microstructure to the global conductance. It is shown that the global 
conductance reduction factor is the product of the tortuosity squared divided by the constriction 
factor, thereby proving that the combined effect of tortuosity and constriction, in addition to the 
porosity and conductance of the electrolyte, fully describes the effective electrical conductance of 
a porous medium. We show that our tortuosity, constriction factor, and conductance reduction 
factor reproduce the electrical conductance for idealized porous media. They are also applied to 
Bentheimer sandstone, where we describe a microstructure-related correlation between porosity and 
conductivity using both the global conductance reduction factor and the distinct contributions from 
tortuosity and constriction. Overall, this work shows how the empirical Archie cementation expo¬ 
nent can be substituted by more descriptive, physical parameters, either by the global conductance 
reduction factor or by tortuosity and constriction. 

PACS numbers: 47.56.-|-r,91.60.Tn 


I. INTRODUCTION 

The electrical conductivity of a brine filled sedimen¬ 
tary rock is an essential part of electrical resistivity log 
interpretation in petroleum exploration [I|. Its impor¬ 
tance has led to a substantial amount of research, with¬ 
out resulting in any satisfactory theoretical description 
of the connection between the pore microstructure and 
the electrical conductivity of the rock . 

The electrical conductivity of a porous medium (e.g., a 
sedimentary rock) consisting of an insulating matrix and 
an electrolyte (e.g., brine) with constant conductivity de¬ 
pends upon the electrical conductivity of the fluid, rela¬ 
tive volumes of the matrix and the fluid (porosity), and 
upon the geometrical distributions of the matrix and the 
fluid (pore microstructure) relative to the applied elec¬ 
trical potential 0. The formation resistivity factor F 
is defined to be the ratio of the electrical resistivity of 
the rock saturated with brine Ro to the resistivity of the 
brine Rw- 

IXw 

In Q Archie formulated an empirical relationship be¬ 
tween the formation resistivity factor F and the porosity 
(j) to describe the conductivity of a porous medium. This 
relationship was given by the equation F = where 

TO is the cementation exponent. 

We will consider the electrical conductivity of porous 
media filled with an electrolyte, or likewise the conduc- 
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tance of networks of resistors. A network might represent 
a simplification of the pore microstructure of a porous 
medium, where pore bodies corresponds to the nodes of 
the network and pore throats to the resistors Q. Net¬ 
work analogues has been widely used to rmresent natural 
porous media such as reservoir rocks [sHtI. 

Several definitions for tortuosity exist. A common def¬ 
inition of the tortuosity r for a path is the ratio between 
the distance I between the endpoints of the path and the 
path length s, i.e. r = l/s (see e.g., Q). It is a measure 
of deviation from the shortest possible path. For an ide¬ 
alized porous media consisting of a single circular tube 
of constant cross-sectional area and length s connecting 
the opposite sides of a cube of side length AS we have 
F = 1/(t^ 0), where (f is the porosity and t — AS/s is 
the tortuosity of the tube 0. 

In an electrolyte the mobile charged particles move 
randomly, however the drift (average) velocity v is given 
by the equation v = fiE for an electro static field E, 
where p, is the electrical mobility. Since E = — V$ for 
the electrical potential $, we have v = —/iVd). For ex¬ 
ample, for a straight circular tube the drift velocity v 
is proportional to / j A, where / is the current and A 
the cross-sectional area. Varying the cross-sectional area 
gives a varying drift velocity, hence energy is expended 
and the conductivity of the tube is reduced. Assume we 
have a tube of length I and cross-sectional area A{x) at 
point X that is filled with an electrolyte with constant 
conductivity a. The conductance of the tube is then 
Gr = <J! f^\IA(x)dx. A tube with equal volume and 
equal length, but with a constant cross-sectional area, 
would have a cross-sectional area A = Jg A{x)dx/l and a 
conductance Go = a A/I = a ^^A{x)dx/P. The inverse 
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of the reduced conductivity due to variation in the cross- 
sectional area is given by the constriction factor [^, [l^ 

C = Go/Gr = 4 / A{x)dx [ -j^dx. (1) 
‘ Jo Jo 

For an idealized porous medium consisting of a single 
straight circular tube with varying cross-sectional area, 
we then have F = C/(f>. 

Since the current is constant in the pore channel, we 
have v{x) oc IjAix), so 

/■' 1 

A{x)dx / , , . dx 

Jo A{x) 

1 

v(x)dx / , , dx 

Jo v{x) 

In the following we use Eq. @ as our definition of the 
constriction factor, which is then a reformulation of the 
more common definition given by Eq. 0 Bn,®- This 
definition also holds when the cross-sectional area is ill 
defined. The constriction factor is a measure of changes 
in drift velocity. 

Note that variation in cross-sectional area also intro¬ 
duces tortuosity. Eor simplicity we disregarded the tor¬ 
tuosity introduced by constrictions when we related the 
constriction factor to the formation factor above. Like¬ 
wise, tortuosity might introduce changes in drift velocity; 
for simplicity this was disregarded when we related tor¬ 
tuosity to the formation factor above. 

The aim of this article is to investigate the physical 
relation between the conductance of a porous medium 
filled with an electrolyte and its microstructure. In Sec. 
m we introduce the concept of the conductance reduc¬ 
tion factor of a porous medium or a network of resistors, 
which is used to describe the local contribution to the 
global conductance. In Sec. IIIII we show how this global 
condnctance can be described by distinct contributions 
of tortuosity and constriction, both inherent to the mi¬ 
crostructure of the porous medium. The relation among 
the global conductance reduction factor, tortuosity, and 
the constriction factor is demonstrated through idealized 
porous media in Sec. IIVI In Sec. |V] we use Bentheimer 
sandstone data to demonstrate these relations also for 
natural porous media. 



II. CONDUCTANCE REDUCTION EACTOR 
DERIVATION 


the conductance reduction factor as 


6 


||V$|| 


As 

a¥' 


( 3 ) 


As mentioned, the potential gradient gives rise to an 
average drift velocity of the mobile charged particles. 
The drift velocity is a conservative vector field, and the 
electric field lines are paths that are everywhere paral¬ 
lel to the drift velocity (equivalently, the electric field). 
Consider an electric field line of length As between in¬ 
let and outlet and with constant drift velocity. Then 
||V<I>|| = A$/As, so 6 = I on the given electric field line. 
In the following we refer to a pore geometry giving a con¬ 
stant conductance reduction factor of t = I as optimal 
(conducting) geometry. A set of straight circular tubes 
oriented in the direction of the applied potential and with 
a constant cross-sectional area gives an optimal geome¬ 
try. An optimal geometry has maximum conductance for 
the given pore volume (porosity), and the conductance is 
^ j Fq — (j) j Rnj j so 


F = 


Rn, 


( 4 ) 


For more complex pore geometries the factor l is no 
longer constant in general. We define the global condue- 
tance reduction factor as the volume weighted average 
of the local conductance reduction factor : 

= (5) 

In the following we show that for any porous medium 
we have F = l/{i^(j>) > l/cf. Hence = (j)^, where 
m is the Archie cementation exponent. Note that the 
conductance reduction factor i might be larger than 1, 
while tg < 1- 

Following an electric field line from inlet to outlet, the 
potential drop along the field line equals the total po¬ 
tential drop over the sample; hence the potential can be 
treated as a measure of the distance between inlet and 
outlet. Comparing the distance on a electric field line 
(ds) for a given potential drop (d$) to the distance in 
an optimal geometry given the same potential drop gives 
the fraction 

dd) As 
ds A$ 

The limit ds ^ 0 then yields Eq. We will see later 
that the conductance reduction factor measures the com¬ 
bined effect of tortuosity and changes in drift velocity on 
the conductivity of the porous medium. 

We next investigate the conductance reduction factor 
for porous media and then turn our attention to network 
analogs. 


Consider a cubed porous medium V of side length 
As consisting of an insulating matrix (without surface 
conductance) and a pore space U C U filled with an 
electrolyte with constant conductance cr, and a potential 
drop A$ over two opposite sides of the cube. We define 


A. Porous media 

For a cubed insulating porous medium of side length 
As hlled with an electrolyte with constant conductivity 
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FIG. 1. An infinitesimal volume element d^l where ril • V"!* = 0 
and nl • Vtl- = iHV-FH. 


a, let a potential drop A$ be applied over two oppo¬ 
site side planes of the cube. For a point x inside the 
pore space assume an infinitesimal volume element 
dfl around x as in Fig. [T] We will orient dfl so that the 
normal vector for the side planes ris is orthogonal to the 
gradient of the potential drop, Us ■ — 0, hence the 

volume is oriented so that there is no current over the 
side planes. Also the normal vector for the end planes rie 
is in the same or opposite direction as the gradient of the 
potential drop, rTe ■ V$ = ±||V<i>||. The potential drop 
between the end planes of dil will be denoted d^. When 
dn —i- 0 the lengths of the electric field lines contained 
in dn converges, and their length will be denoted ds. 
Also the cross-sectional areas in dil will converge when 
dft —>■ 0, and will be denoted by dA. The resistance over 
dQ is then given by dR = ds/{adA). 

We will now investigate the corresponding volume dD! 
in an optimal geometry that gives equal current dl' = dl 
given the same potential drop d$' = d$, where all 
marked parameters are from the optimal geometry. In¬ 
voking Ohm’s law gives dR' = d^'/dl' = d^/dl = dR. 
Since t' = 1 in the optimal geometry, we have ds' = 
Asdd>'/A<I) = Asd$/A<i). Using 


ds' 

adA' 


= dR' = dR = 


ds 

adA ’ 


we then obtain 


medium with pore space H. Invoking Eq. (j4]) this yields 
V V 1 

n' ilf ^ ’ 

where V is the total volume of the porous medium. Hence 
we have derived a relationship between the local conduc¬ 
tance reduction factor and the global conductance reduc¬ 
tion factor that describes the global conductance. 

We might assume that for every point x € we have 
continuous paths in Q to the two sides where the po¬ 
tential drop is applied. All parts of the pore space not 
satisfying this condition can be assumed to be part of the 
matrix, since it will not contribute to the conductance. 
However, since i^{x) = 0 for all points x € not satis¬ 
fying the condition above, the relation F = l/{i^(j)) still 
holds when such points are included in H. 


B. Network model 


Consider a network representation of a porous medium 
consisting of resistors (pore throats) and nodes (pore 
bodies) connecting the resistors A resistor t has an 
associated length s{t) and a constant cross section A{t). 
The nodes do not have associated volumes, they are only 
connecting the resistors. In such a network consider a 
resistor t of volume V{t) = A{t)s{t), and with potential 
drop over its length A$(t). The resistance over t is then 
R{t) = s{t)/{aA{t)). For the whole sample we denote 
the potential drop A<i> and the size of the sample in the 
direction of the global potential drop As. 

We now investigate the corresponding volume V'{t) in 
an optimal geometry that has equal current A$(f)/i?(t) 
given the same potential drop A$(t). Then R'{t) = R{t), 
which leads to 


s'jt) 

aA'{t) 


R'{t) = R(t) = 


s{t) 

aA{t) 


dA'ds' 


ds'dA , dA 
ds ® ^ 
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dAds. 


Since dA'ds' = dU', we obtain 


,, dU' / dd) As 

lim —— = lim ——-— 
dQ —^0 Clol dri —^0 V ds 


= l|V<i>(x)|| 


AsV 


= i 


This leads to 


n' = [ dn' = [ t^du = tin, (6) 

Jq Jn 

where is the volume needed in an optimal geometry 
to obtain the same electrical conductance as the porous 


Using that 6' = 1 in the optimal geometry, we have 
s'(t) = A$(t)As/A$. Then, similarly to the porous 
medium, we find 


A'it)s'{t) 


s{t) 


At) 


m 

s{t) 



2 


V s{t) 



A{t)s{t) 


A{t)s{t). 


Hence V'{t) = A'{t)s'{t) = A{t)s{t) = t(<)^U(t) is 
the volume needed in optimal geometry to obtain equal 
current given the same potential drop as in resistor t. 

Let n — EU(t) and = T,V'{t) be the total network 
volume and the volume needed in an optimal geometry to 
transport the same current, respectively. For a network 
we define 
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Hence the global conductance reduction factor for the 
whole sample is a volume weighted average of the local 
conductance reduction factors With il' = and 

invoking Eq. we have 

V V 1 

F = _=_=_ ffll 

n' ijn Lif ^ ’ 

where V is the total volume of the porous medium. 

It is important to note that the local components of 
the global conductance reduction factor are inter-related 
since the local potential drop is dependent on the full 
network. Nevertheless, the conductance reduction factor 
offers the possibility to investigate the local contributions 
to the global conductance. 


III. TORTUOSITY AND CONSTRICTION 


It is well established that tortuosity and constriction 
affect the conductance of a porous medium 0, i[ni- 
What we show next is how tortuosity and constriction 
affect conductance. 

As before, consider a cubed insulating porous medium 
of side length As filled with an electrolyte with constant 
conductance and an applied potential drop A$ over two 
opposite sides of the cube. We define the tortuosity t of 
this porous medium as the side length of the cube divided 
by the average electric field line length. This is given by 


r = 




1 

As 


( 10 ) 


where H is the pore space and dU is an infinitesimal vol¬ 
ume element oriented in alignment with the drift velocity 
(see Fig. [T]), dl is the current transported through dO, 
and It is the total current through the porous medium. 
Here we use that 


lim 

dn^o 


dids 

dU 


lim = c^l|V«>||. 

dn-^o ds 


Note that ds is the length of the infinitesimal dft in the 
direction of the drift velocity, so the integral in Eq. m 
yields a current-weighted average of lengths in the direc¬ 
tion of the drift velocity. In other words, the integral 
gives the average distance the current has to traverse. 

Since we introduce a measure for constriction along 
electric field lines, we need a measure for tortuosity along 
electric field lines too. Let E € P be an electric field line, 
where P is the set of all electric field lines, let Ir be the 
length of the electric field line E, and let He = {x G 
Hlx € r € P} = {x G H|V<I>(x) ^ 0} be the subset of 
H given by the set of electric field lines or, equivalently, 
the conducting subset where V<I> 7 ^ 0. We define the 
conducting porosity as the fraction (pc = Ilc/V. The 
tortuosity Tc of the porous medium is given by 


1 

He 


r(r) 2 dr = ^ / 



( 11 ) 


Here T(r) = As/Zr is the tortuosity of the electric field 
line T. We can, without loss of generality, assume that 
the set dr is a simply connected space (i.e., there are no 
holes in dT). This assumption is used throughout this 
article. 

The constriction factor was introduced in Sec. U and 
measures changes in drift velocity. The constriction fac¬ 
tor C'(r) for an electric field line T of length Ir is given 
by the integral 




ds 




1 


||V$(s)|| 


ds. 


( 12 ) 


where we use that lim^is^o = ||V<i)|| on the elec¬ 

tric field line T. We define the constriction factor Cc 
as the current-weighted average of the electric field line 
constriction factors C'(r): 


Cc = Y [ C{T)dIr. (13) 

h Jv 


Here d/r is the infinitesimal current to the infinitesimal 
set of electric filed lines dT, and as before It is the total 
current through the porous medium. 

We then obtain Lc by restricting ig as given by Eq. m 
to lie- 




1 


H,, 


/ t^dH = 


t(r)2dr 




(14) 


Here ^(r)^ = lim^r-^r l/<^r is the conductance 

reduction factor of the electric field line F. Note that 


L 


2 

c 




(15) 


since i{x) = 0 for x G H \ He. 

By choosing all dH such that they contain all electric 
field lines in dT, we have the following correspondence: 


t(r)^ = lim 


1 


dr^r dr 


i^dH 


= lim 


'dr 

1 


d<I>\ /As 


dr^r JpdA(s)ds jY-yds) \A$ 
As \ ^ 1 f f d^\^ ds 


dA{s)ds 


ds 


IrmdsJr\dsJ d4> 

T(r)= 




As 


c(r). 


(16) 


Here we use the correspondence dip = adA{s)d^/ds, 
where as before d/p is the infinitesimal current to the 
infinitesimal set of electric field lines dT. Also, as before 
the set dr is simply connected. The cross-sectional area 
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of dT at point s along the electric field line T is denoted 
dA{s). 

We now have the framework needed for our main re¬ 
sult; combining Eqs. (HJ), (HSl), and m we obtain 


ft tfrCmdir “/,C(r).(r)2<ir “■ ^ ' 

Here we employ the equality 

^ r2cAHA$(.^(T 

‘ ^ A^ ^ EA^ 

and the corresponding equality 

drAHA$6(r)2cr 

dJr = -:- 

EAs. 

For a porous medium, combining Eqs ©, (USD, and 
(113, we then arrive at the equality 


F=± = ^ = ^. 

T^4>c 


(18) 


The constriction factor Cc measures the loss of conduc¬ 
tance due to changes in drift velocity, and the tortuosity 
measures the loss of conductance due to electric field 
line lengths. Given by the equality above, these two fac¬ 
tors fully describe the electrical conductance for a given 
porosity (pc and conductance of the electrolyte a. The 
tortuosity and constriction factor are describing the pore 
microstructure relative to the applied electrical potential, 
thus they give a theoretical description of the connection 
between the pore microstructure and the electrical con¬ 
ductivity of the porous medium. 

If we assume that for every point a; G we have con¬ 
tinuous paths in H to the two sides where the potential 
drop is applied, then V<i>(a;) ^ 0 almost everywhere in 
n, so the closure of the set ilc equals fl. However, so¬ 
lutions for the electric potential in a grid representation 
of the porous medium might still give V<I> = 0. Also, 
in resistor network analogs where every element has net¬ 
work paths both to the inlet and to the outlet, we might 
have V$ = 0, for example, when nodes are connected to 
exactly one resistor. Hence flc C fl in general. However, 
when (p = (pci we have 


1 




(19) 


We perform calculations of tortuosity Tc and constriction 
factor Cc on resistor network analogs of porous media in 

Sec. El 



FIG. 2. An idealized porous medium consisting of one tube 
of length s with constant cross-sectional area inside a cube of 
side-length AS. 


AS, and with a potential drop A$ over the side planes 
connected by the pore channel. For all examples in this 
section we then have Uc = ^i (pc = (pi = t-gi Cc = C 
and Tc = T. Therefore we do not use the c subscripts. 


A. Conductance reduction factor and tortuosity 


Consider a tube of length s and constant cross- 
sectional area, as depicted in Fig. [5] The magnitude 
of the gradient of the potential inside the pore channel is 
constant at A$/s, thus 




Since b is constant throughout the pore channel, we ob¬ 
tain 


1 _ 1 _ 1 
Lgp L^P t'^P 


( 20 ) 


B. Conductance reduction factor and constriction 

In this subsection, we consider an idealized model of 
constriction. Consider one straight circular pore channel 
in the direction of the potential drop, hence of length 
AS, with cross-sectional area A(s) for s G [0, AS']. In the 
following we neglect the tortuosity introduced in such a 
medium due to the constriction, as mentioned in Section 
m Since the electrical current I = I{s) = aA{s)^ = 
crA(s)||V$(s)|| is constant, we get 


IV. EXAMPLES FROM IDEALIZED POROUS 
MEDIA 


rAS 


rAS 


A$ = 


d^ = — 


ds 

W)' 


Let us return to the idealized porous media described 
in Sec. HI consisting of one circular pore channel con¬ 
necting to opposite side planes of a cube of side-length 


||V$(s)|| 


I 

aA{s) 


A^ 

Ms) l/A{s)ds 





















6 


Then the conductance reduction factor is given by 

= ||V4.(s)f 

_y. 

\Ms) ^/Ms)ds J 

By volume averaging the conductance reduction factor 
i(s)^, we have 


rAS 




A{s)ds 


rAS 


AS'^ 


C'A{s)dsJo (Ais)) 

AS^ 


-A{s)ds 


Jo ^ A(s)ds /p ^ds 


= 1/C. 


This yields 



C 


( 21 ) 



FIG. 3. A visualization of the network representation of the 
pore space of Bentheimer rock model / with 15.4% porosity. 
Balls represent pore bodies, and sticks represent pore throats. 
In this network there are 9260 pore bodies and 14792 pore 
throats. 


C. Conductance reduction factor, tortuosity, and 
constriction 


Let us turn to a combination of tortuosity and con¬ 
striction. We now consider a single circular pore channel 
with cross-sectional area A(s) for s G [0,s], where s is 
the length of the tube. In line with the calculations in 
the previous subsection, substituting the length of the 
channel with the new length s, we obtain 


c(sf = 


AS 


i/A{s)ds 


^1 = 


AS^ 


J^Ais)dsJo ^dsr(A(s)y 


A(s)ds 


AS 


thus yielding 


s J fj A(s)ds ^ds 


= r^C, 


F = 


1 C 

Fg(l) 


( 22 ) 


V. BENTHEIMER ROCK EXAMPLE 

In this section we use a micro-CT (microtomography) 
image and rock models of a Bentheimer sandstone to 
compute the global conductance reduction factor t^, the 
tortuosity Tc, and the constriction factor Cc- Using the 


commercially available software e-Core we generated 
three-dimensional rock models of a Bentheimer sandstone 
with porosities ranging from 5% to 25%. The software 
simulates the processes of grain sedimentation and dia¬ 
genesis [l^ . The same grain sedimentation model pa¬ 
rameters were used for all models. We achieved different 
porosities by changing the amount of quartz cementation 
|I4l| . as occurring in diagenetic rock forming processes, 
while other processes like compaction could also natu¬ 
rally lead to a porosity reduction. The model sample 
size was 2.5 mm^ with a resolution of 5 ^m, compared to 
the micro-CT image of 3.7x3.7x2.4 mm with a resolution 
of 3.7 /rm. 

Network representations of the pore microstructure 
were numerically extracted from both the generated mod¬ 
els and the micro-CT image. The algorithm segments the 
grid representation of the pore microstructure into pore 
bodies and pore throats, each associated with a volume 
and a point x G V. In Fig. [3] we have visualized one 
such network. The distance ||a;ti — Xt 2 \\ between the two 
pore bodies tl and t2 connected by a pore throat t3 is 
divided into parts: parts In and lt 2 associated with each 
pore body, tl and t2, respectively; and part Its associated 
with the pore throat, where 1 ( 14 -/ 43 -|-Zt 2 = ||a^ti—a:t 2 ||- We 
associate tubes with a constant cross-sectional area with 
these parts: For the middle tube the cross-sectional area 
is Vts/ks with length Its, while for the tubes associated 
with the pore bodies tl and t2 the cross-sectional areas 
are Vu/(atihi) with lengths lu, where Vu is the volume of 
the pore body or throat as given by the segmentation al¬ 
gorithm, and ati is the number of pore throats connected 
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to pore body ti. These tubes then have conductances 
Git = (^Vtilil'^^ati) for i = tl,t2 and Gts = aVajlt^. 

We next transform the network model to a resistor 
network analog. There is a one-to-one correspondence 
between the pore throats in the porous medium and the 
resistors in the resistor network analog; also, there is a 
one-to-one correspondence between the pore bodies and 
the network nodes. Each resistor t is given a conductance 


Gt = 



fl]^ ^ l]^\ 

\ Vtl ^ Vt 2 J ’ 


(23) 


where subscripts t\ and t2 represent the pore bodies con¬ 
nected to throat t3. In contrast with e.g., Eq. (18) in 
Ref. [l^, this conductance formulation ensures that the 
resistor network analog has the same pore volume as the 
corresponding rock model or micro-CT data. 

Let be the potential corresponding to node and 
tij the resistors connected to node i. We then solve for 
<i>i such that 


(Jii 


t=i 


$,) = 0, 


(24) 


where inlet and outlet nodes have fixed potential as 
boundary conditions. This is equivalent to Kirchhoff’s 
circuit laws for the conservation of charge and energy in 
electrical circuits. The solution technique is equivalent 
to solving a resistor network problem. 

The resistor t is viewed as consisting of the three 
tubes and t2 with a constant cross-sectional area 

in series. We find potentials for i = 1,2 such that 
— <l>it| = Gt 3 |<I>it — $ 2 t| = Gt 2 \^ 2 t — ^t 2 \- This 
enables the calculation of the conductance reduction fac¬ 
tor for the three tube sections: 



The volume average of these local contributions then 
gives the global conductance reduction factor. 


1 




t 




tl 2 


Vt 


4 + —4 + 


Oil 


0*2 


(25) 


The calculated results are reported in Table HI 

In Table U we have also listed the formation resistivity 
factor as F = l/(tg())). Fig. |4] is a plot of the porosity 
and formation resistivity factor for the network represen¬ 
tations. Here the micro-CT data fall in line with the rock 
models. 

In Fig. [5] we have plotted the global conductance re¬ 
duction factor ig versus the porosity (j). The linear func¬ 
tion = 1.010(())— 0.018) is included as a trend line 
for the data. For (j) = 0.018 we have = 0, which is 



Porosity (fraction) 


FIG. 4. Plot showing porosity tj) versus formation resistivity 
factor F for both rock models and micro-CT (mCT) data, 
together with two functions describing the correlation. 



FIG. 5. Plot showing global conductance reduction factor 
versus porosity <j) for the rock models and the micro-CT 
(mCT) data, together with a linear fit to the data. 


interpreted to be the percolation threshold for this sand¬ 
stone [HI, [l^. Thus 

tg((/))(/) (/>(</) — 0.018) (j>{4) — 0.018) 


This function is also included in Fig. 0] The match with 
the data is fair, given the simplicity of the function. 

The network volume with nonzero potential gradient 
can be calculated as 


He - E 



Ott2 




t3 


(27) 


In Fig. [n]we have plotted porosity 4> versus conducting 
















TABLE I. Model and micro-CT results. 



Porosity 

Gonduction reduction factor 

Tortuosity 

Constriction 

Formation resistivity 

Gementation 

Rock 


0C 



T 

Tc 

Cc 

factor F 

exponent m 

a 

0.049 

0.028 

0.017 

0.030 

0.430 

0.399 

5.290 

1186.03 

2.355 

b 

0.076 

0.053 

0.038 

0.054 

0.471 

0.441 

3.574 

347.18 

2.272 

c 

0.096 

0.077 

0.083 

0.103 

0.542 

0.519 

2.604 

125.80 

2.064 

d 

0.122 

0.106 

0.116 

0.133 

0.569 

0.546 

2.241 

71.02 

2.024 

e 

0.138 

0.124 

0.135 

0.151 

0.588 

0.567 

2.126 

53.43 

2.011 

f 

0.154 

0.139 

0.150 

0.165 

0.595 

0.576 

2.011 

43.51 

2.014 

g 

0.159 

0.146 

0.153 

0.167 

0.602 

0.581 

2.028 

41.03 

2.020 

h 

0.177 

0.166 

0.172 

0.184 

0.609 

0.589 

1.888 

32.83 

2.016 

i 

0.199 

0.188 

0.193 

0.204 

0.632 

0.611 

1.832 

26.04 

2.020 

j 

0.221 

0.211 

0.203 

0.212 

0.633 

0.612 

1.769 

22.33 

2.056 

k 

0.241 

0.233 

0.217 

0.225 

0.645 

0.624 

1.731 

19.07 

2.073 

1 

0.254 

0.247 

0.219 

0.226 

0.646 

0.624 

1.726 

17.97 

2.109 

Micro-GT 

0.193 

0.169 

0.171 

0.195 

0.629 

0.612 

1.924 

30.40 

2.075 



Porosity (fraction) 

FIG. 6. Plot showing the correspondence between the poros¬ 
ity (j) and the conducting porosity cpc. for the rock models and 
the micro-CT (mCT) data, together with a linear fit to the 
data. 

porosity (j)c = ^c/V. A linear fit to the data in Fig. [6] 
gives the correspondence 

^^(<^) = 1.07(/>- 0.025, (28) 

which we use in the following. For </> = 0.025/1.07 = 
0.023 we have 4'c{4>) = 0, which is interpreted as a perco¬ 
lation threshold. This is compared to the threshold 0.018 
derived from the global condnctance reduction factor in 

Eq. dUD. 

For each model and the micro-CT image we discretized 
the volume flc into a disjoint union UF, where volume F 
conducts a current Ip such that '^It = It- Each volume 
F is contained in a single series of resistors, with the 
first resistor connected to an inlet boundary and the last 
connected to an outlet boundary. Since we cannot trace 
current across the network nodes in a porous medium, 



(T ortuosity)^/Constriction 

FIG. 7. Plot showing the correspondence between the con¬ 
ductance reduction factor t(r)^ and the tortuosity squared 
divided by the constriction r(r)^/C(r). Each dot represents 
a volume F in the discretization flc = UP for the resistor 
network analog of the micro-GT image of a Bentheimer sand¬ 
stone. 


the discretization flc = LIF is not unique. With respect 
to a resistor network analog, full mixing is implied at the 
node points. Different algorithms for the discretization 
were tested on smaller network analogs for porous media, 
all giving results comparable to the discretization derived 
from full mixing at the nodes. 

The parts {Fti} of F is in tubes i = 1,2,3 of con¬ 
stant cross-sectional area as associated with resistor t, 
as described before. Each volume F^i C Vu then has 
the associated length la- For each F we calculated the 
conductance reduction factor 

= p X/ + ^‘ 2^42 + ^43643), ( 29 ) 

4 
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FIG. 8. Plot showing porosity (p versus both tortuosity 
squared and inverse constriction for the rock models 
and the micro-CT (mCT) data. 


the tortuosity 


r(r) = 


As 


+ ^42 + ^43 ’ 

and the constriction factor 
A$ 


c(r) = 


(S 4 + ^*2 + ha) 


E 


/2 

hi 


l$41 - $1 


+ 


^t2 


1^/2 — ^2 


(30) 


43 


|$1, _ $2 


(31) 


From Fig. [7]we see that (-(F)^ equals r(F)^/(7(r) for all 
F, hence the values from Eqs. (EH)) to (EH) are consistent 
with Eq. (HU). 

We then obtain the tortuosity and constriction factor 
for the volume ^c'- 



(32) 

ft = |y;/rC(r). 

— 

(33) 


The calculated values are listed in Table ID We see that 
T^/Cc = which is consistent with Eq. (fT71) . The 
partition of the global conductance reduction factor into 
tortuosity and constriction gives additional information 
on the electrical conductance compared to the traditional 
Archie theory. Note that the values for Tc and r are 
similar. From Table |T] we observe that (f>i^ = (pc^c = 
(pcTc jCc, which is in correspondence with the theoretical 
derivation in Eq. EH). 

In Fig. [5]we have plotted the porosity p versus both the 
tortuosity and the inverse constriction C”, together 
with functions fitted to the data. For high porosities, i.e., 


rocks with little cementation, we have less conductivity 
reduction due to both tortuosity and constriction, com¬ 
pared to that for low porosities. This is as anticipated, 
as cementation increases the ratio between the pore body 
and the pore throat cross-sectional area, which yields a 
larger fluctuation in drift velocity and hence a higher con¬ 
striction factor. Cementation also can block pore throats 
completely, which increases the length of the electric field 
lines and hence reduces the tortuosity. With increasing 
amounts of cementation (decreasing porosity), tortuos¬ 
ity becomes as important as constriction in reducing the 
conductivity. 

The tortuosity of Bentheimer sandstone is matched 
fairly well with the function t^((/>) = 0.150ln(())) -|- 
0.609, and the constriction with the function C~{4>) = 
0.244 ln((^) -I- 0.937, both plotted in Fig. [5] The func¬ 
tion for the constriction indicates a percolation threshold 
of = 0.020, while the tortuosity indicates a 

percolation threshold of e-0-609/o.i50 _ q Qlg, which is in 
line with the result derived from the global conductance 
reduction factor in Eq. (1^51) and with the result from the 
conducting porosity in Eq. (1^51) . Combining Eqs. EH 
and (l33l) for tortuosity and constriction with Eq. (l28l) 
yields 


m = 


'I{4>)Cc {4>)4'c 


1 


(0.15 ln(()>) -b0.61)(0.241n(()i) -b0.94)(1.07(/>- 0.025) 

26-3 

(ln(((.) -b 4.02)(ln((^) -b 3.93)(()) - 0.024) ^ ^ 


This relation is also plotted in Fig. |4] and provides a 
better match to the Bentheimer data than Eq. (EH). 


VI. CONCLUSION 

In this work we have described and calculated the elec¬ 
trical conductance of an insulating porous medium filled 
with an electrolyte. We have shown how the conduc¬ 
tance depends on the length of the electric field lines and 
changes in drift velocity, in addition to the porosity and 
the conductance of the electrolyte. The lengths of the 
electric field lines are described by the tortuosity Tc, and 
changes in drift velocity are described by the constriction 
factor Cc, both inherent to the microstructure of a porous 
medium. The combined effect of the tortuosity and con¬ 
striction is defined as the global conductance reduction 
factor When pc = P, we have 6^ = = r^/Cc, where 

Lg is the volume weighted average of 



We have shown that our methodology reproduces re¬ 
sults for idealized porous media. Our derivation is 
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also applied to natural porous media given by Ben- 
theimer sandstone, where we describe a correlation be¬ 
tween porosity and conductivity using the global con¬ 
ductance reduction factor and the distinct contributions 
calculated for tortuosity Tc and constriction Cc- 

This work enables the substitution of the Archie ce¬ 
mentation exponent with a more descriptive, physical 
parametrization; either with the global conductance re¬ 
duction factor or with the tortuosity and constriction fac¬ 
tors, separately. The latter parameters rely on accurate 


three-dimensional information on the porous medium mi¬ 
crostructure, e.g., micro-CT images. Fortunately, such 
microstructure data become increasingly more available. 
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